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Abstract 

We include the new, five-loop, O(a^) correction into the QCD sum rule used 
for the s-quark mass determination. The pseudoscalar Borel sum rule is taken as a 
study case. The OPE for the correlation function with N^LO, O(af) accuracy in the 
perturbative part, and with dimension d < 6 operators reveals a good convergence. 
We observe a significant improvement of stability of the sum rule with respect to 
the variation of the renormalization scale after including the O(af) correction. We 
obtain the interval m s {2 GeV) = 105 ± 6 ± 7 MeV, which exhibits about 2 MeV 
increase of the central value, if the O(a^) terms are removed. 
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1 Introduction 


A precise determination of the strange quark mass m s is extremely important for various 
tests of Standard Model. A reach variety of approaches is used to evaluate this funda¬ 
mental parameter in QCD. Recently, the first unquenched lattice QCD determinations 
became available [1-6]. In addition, ChPT provides rather accurate ratios of strange and 
nonstrange quark masses [7,8]. Furthermore, one evaluates m s , combining the operator 
product expansion (OPE) of various correlation functions for strangeness-changing quark 
currents with dispersion relations. These methods include model-independent bounds [9], 
QCD analyses of hadronic r decays (see [10-14] for the latest results), as well as different 
versions of QCD sum rules [15] and related finite-energy sum rules (FESR) [16]. The most 
recent sum rule determinations of m s in the channels of scalar, pseudoscalar and vector 
currents are presented in [17-19], respectively, references to earlier analyses can be found 
in reviews [20,21] (see also [22]). The estimated accuracy of these results is 15 -30%. To 
achieve a better precision, one has to calculate higher orders in OPE of the correlation 
functions and gain a better control over potentially important nonperturbative corrections 
beyond OPE (the so called ” direct instantons”). Furthermore, more accurate data for the 
inputs in hadronic spectral functions and a better assessment of the quark-hadron duality 
are needed. 

In this paper we concentrate on the QCD sum rules used to evaluate the strange quark 
mass and make one further step to improve the accuracy of this determination by including 
the N 4 LO perturbative QCD corrections of 0(aj) into the sum rule. The 0(a 4 ), Eve-loop 
contribution has recently been calculated for the correlator of the scalar quark currents in 
[23] and can be equally well used for both scalar and pseudoscalar sum rules. As a study 
case, we choose the pseudoscalar version of the standard Borel sum rule. For the hadronic 
spectral function we employ the three-resonance ansatz worked out in [18]. 

We find that both OPE for the correlation function and the resulting sum rule reveal 
a good numerical convergence in powers of a s . The new 0(a 4 ) correction to the sum rule 
has a naturally small influence, resulting in about 2 MeV decrease of the s-quark mass m s 
(in MS scheme) determined with O(a^) accuracy. Importantly, after including the O(af) 
correction, we observe a significant improvement in the stability of the extracted value 
of m s with respect to the renormalization scale variation in the sum rule. With 0(a 4 ) 
accuracy we obtain the interval 


m s (2 GeV) 


105 ±6 


±7 


param 


hadr 


MeV, 


( 1 ) 


where the estimated uncertainties from the sum rule parameters and hadronic inputs are 
shown separately and will be explained below. 

In what follows, after a brief recapitulation of the pseudoscalar Borel sum rule in 
Sect. 2, we present in Sect. 3 the QCD OPE expressions for the underlying correlation 
function, including the new 0(a 4 ) terms in the perturbative part. The Borel transform 
and the imaginary part of the correlation function are also given. In Sect. 4 we turn to 
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the numerical analysis of the sum rule and obtain the interval (1). Sect. 5 contains the 
concluding discussion. 


2 Pseudoscalar Sum Rule 


We consider the correlation function: 

n (5 V) =ijd 4 x e iq ' x ( 0 \T {#(z)# f (0)} |0> ( 2 ) 


of the two pseudoscalar strangeness-changing quark currents, defined as the divergences of 
the corresponding axial-vector currents: 

j = &*(s7t*75q) = (m s + m q )si^ 5 q. (3) 


For definiteness, the light quark q = u is taken. 

The Borel sum rule is obtained following the standard SVZ method [15] and is based on 
the (double-subtracted) dispersion relation for Il (5 ^(g 2 ). This relation is more conveniently 
written in a form of the second derivative: 


n ( 5 ) V) - 


d{q 


After Borel transformation one obtains: 2 


2 
7 T 



i m n ( 5 )(s) 
(s — q 2 ) 3 


n (5 ) "(m 2 ) 


Bm 2 


n (5) "(g 2 ) 


/ *e“* /Ml Imn( 5 >(s). 
0 


( 4 ) 


( 5 ) 


The l.h.s. of the above relation is calculated in QCD at large M 2 S> A q CD i n a form 
of OPE (perturbative and condensate expansion), in powers of a s and m s /M , and up 
to a certain dimension of vacuum condensates. In r.h.s. the hadronic spectral density 
PhldM = (l/7r)Imrk 5 )(s) is substituted. At s < so , where so is some effective threshold, 
the function p^l dr ( s ) includes kaon and its excitations. The rest of the hadronic dispersion 
integral at s > so is approximated using quark-hadron duality, p^l dr {s) ~ Pope ( s ), with 
the spectral function calculated from OPE: Pqpe ( s ) = (1/7r)Ini [n (r 4(s)]op.E- The final 
form of the sum rule is: 


SO OO 

M 4 [n« 5 >"(M 2 )]oPE = Jdse-W V® r (s) + Jdse-/ M VghW . (6) 

0 so 

In the following we discuss both parts of this equation in detail. 

2 Here we use the following normalization convention: f3jv/ 2 [l/(a — q 2 )] = e“°/ M . 
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3 OPE results to 0(cfy 

In this section we present the expressions for [II (q 2 )\oPE and, correspondingly, for 
[TlW(M 2 )\ ope and Pqpe( s ) determining the QCD input in the sum rule ( 6 ). The OPE 
for (q 2 )\oPE goes over powers of (1 /q 2 ) d+2 ordered by the dimension d = 0,2,4, 6 . 
The OPE terms with d > 6 are neglected, while already the d — 6 contribution is very 
small in the working region of the variables Q 2 and M 2 . 

The d — 0,2 terms of OPE originate from the perturbative part of the correlation 
function. The expansion in quark-gluon coupling up to four loops, that is, up to 0(a 3 ), 
can be taken from [24-26]. The new O(a^) terms are obtained in [23]. Putting them 
together, we obtain: 


rrr(5) , T^ ) 2u0=0,2) _ 3 {m s + m u 

1 [Q )loPE ~ 8n 2 Q 2 


1 + 


y , «0 ,i 1 



(7) 


where Q 2 = — q 2 , and the coefficients multiplying the powers of the quark-gluon coupling 
a s = a s (p)/n are 


- _ U - _ 5071 35 139 17 2 

“°.i - y - 2 l Q, «o,2 - ~ y & - -g“ l Q + -J l Q, 


( 8 ) 



1995097 1 

5184 36^ 


65869 . 715 . 2720 , 475 _ , 695 , 2 221 ,, , . 

-C 3 3-C 5 -oH- C 3 lo-\ - l n - l n , (9) 

216 12 9 Q 4 y 8 Q 24 Q K J 


do, 4 — 


2361295759 

497664 

52255 

^5^ C? + 1q 


2915 _ 4 25214831 

10368 71 5184 ' 216 

43647875 1 4 864685 

+ XX ^ H XXX S3 


/2 
L Q 


1778273 

1152 


10368 
16785 
32 


18 


192155 2 59875 625 

13 T o-i n ^>3 108 8 "* 43 

24025 
48 


-t -/ 3 

79333' 

1 / 4 

'7735' 


288 . 

+ i Q 

. 384 . 


( 10 ) 


d 2 , 


1 — 



<^ 2,2 


8557 77 147 25 2 

~ti 3 ~ y Q + y Q 


<'u> 


including the new result for do,4. Here Iq = log^-, and ( n = (( n ) is the Riemann’s Zeta- 
function. The coupling a s and the quark masses m s and m u are all taken in MS scheme 
at the renormalization scale /i. We have neglected the light-quark mass m u , except in 
the overall factors. Note also that in the subleading d = 2, O(m^) terms of the above 
expansion, the currently achieved 0(a 2 ) accuracy is quite sufficient. 
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The contributions with d — 4, 6 in the correlation function originate both from nonper- 
turbative (condensate) terms and from 0(m®) corrections, and we use the known expres¬ 
sions [26, 27] 


where 


and 


[n ( 5 >V)l 


2 \i (d=4,6) _ ( m s + m uY 


OPE 


Q 6 


( 23 

-2 m s {uu) I 1 + a s (y - 2 l Q ) 


■i/ c (i + a,(^L - 2 l Q )\ + I. (l + «,(y - 2fqO 

3 4 / 1 155 15, \ h \ 

w m *U + ^“T 1q ) + q->I' 


la — - 


T /— \ 3 4 ( 1 53 

I,=m,( SS ) + — 2 ra, --- 


9 . cx s „ 2 \ (s 16 \ / 91 \ . , 3 / 4 

4<— G ) ( 1 + y a s ) + Aa s ( 1 + ) m s( ss ) + 4^ ( 1 + 3 a " ) rn s 


( 12 ) 

(13) 


(14) 


are the vacuum expectation values of two RG-invariant combinations of dimension 4 con¬ 
taining quark and gluon condensate densities (for details and explanation see [26-29]). 
Finally, 

Iq = —3m s (uuG) ——7 T 2 a s [^(uu) 2 + ( ss ) 2 — 9(uu)(ss)j (15) 

is the combination of dimension-6 contributions of the quark-gluon and 4-quark condensates 
(the vacuum saturation is assumed for the latter). 

The Borel transform of Eqs. (7) and (12) is given by 




1 + b,u a * - 2 jy ( 1 + h 




where Im = log jjr and the coefficients are 
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Vi —7T + 2 7 b — 2 Im, 


r 5071 139 17 2 17 2 35 . 139 , 17 , 17 l2 

J °’ 2 -J 44 + ~ 7E + T 7 E ~Y 7r ~Y Cs ~~ M ~ Y 7E m + T M ’ 
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7e + Ya 7e 


695 


5184 9 
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216 4~ 7e YY 


48 


- 7 r 


221 

48 


7b7T 


+ l 


M 
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' b M 


2720 695 221 

9 t 7e ~yy 

695 221 

■7T" + -ztIe 
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221 

Ya 1m 
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(16) 


(17) 

(18) 


(19) 
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2361295759 | 43647875 _ | 1778273 _ 2 | 79333 3 | 7735 _ 4 1778273_ 2 
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18 288 16 ' 48 J 
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1 ^ ,4 [7735] 

lE \ + 1m [ Wj ’ 


&2,1 — — +47 e ~ 4 Im, 


— 5065 97 25 2 25 2 77 A 97 , _ , 25 , 2 

b-2,2 - + — lE + — - y C 3 - + — l M i (22) 

and, respectively, 

[n (5 >"(M 2 )]gjJ e) = [-2m,{g«> (l + a,(y + 2 7e - 2 <„)] 

— (l + — 3 7 e — 2 /m)] — 2 , [ 1 + a s("g" + 2 7 e — 2 /^)] 

3 _4 /" 1 , 5 , i5 15, ^ /(, 1 

~^ m 'U 6 T 7E “T iM j + 3M7' (23) 

In addition, we need the imaginary part of the correlation function calculated with the 
same aj accuracy as Eqs. (16) and (23): 

fiSlM = hmn® W = 3(m - 8 + m " )2 s 1 1 + £f w <- 2 ^ (i + Y ,*>.<<) } 

, m '( S ) f 45 ^.4/5 , o„ , a ^ s ) T „ 


where l s = log ^ and 


~^2 m s( s ) + 2a a (s)m s (w«) 4 - —I G - a 8 (s)I 8 > , 


17 , 9631 17 2 35 95 , 17, 2 

r °' 1_ T _ r °’ 2 — ~144 12^ _ T C3 ~ Y * + T 
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r 0,3 = 


4748953 


ro,4 = 


229 , 1 4 

- 71 -7T — 

6 36 


5184 6 36 216 

475 A , 229 9 221 , 
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C 3 


(27) 


16 , , _ 5065 

r 2,i ——4t s , r 2,2 — 


25 2 77 . 97 , 25 , 2 

—7T- C3 - h H-C 
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(28) 


The OPE expressions are valid at sufficiently large Q 2 3> A qcd or ) correspondingly, 
at large M 2 . It is well known that in the spin zero (scalar and pseudoscalar) channels 
the breakdown of OPE is expected to occur at relatively large Q 2 ~ 1 GeV 2 , due to 
the presence of nonperturbative vacuum effects which are beyond the local condensate 
expansion [30,31]. Models of the correlation function based on instanton ensembles, such 
as the instanton liquid model (ILM) [32,33] allow to penetrate to smaller Q 2 . A remedy 
used in previous analyses of pseudoscalar sum rules is to add to the OPE series an instanton 
correction calculated in ILM. As realized, e.g., in [18], at sufficiently large M 2 , practically 
already at M 2 > 2 GeV 2 the ILM correction is small, hence we will avoid it by choosing 2 
GeV 2 as a lower limit of the Borcl mass. 

For the reader’s convenience, the lengthy coefficients appearing in eqs. (7), (16) and 
(24) are made available (in computer-readable form) in [34], 


4 Hadronic spectral density and the sum rule 

The spectral function f>^l dr (s) in Eq. ( 6 ) is a positive definite sum of all hadronic states 
with strangeness and J p = 0~, located below the threshold yGo, above which p^dr( s ) i s 
approximated by the OPE spectral density. Clearly, the larger is So, the smaller is the 
sensitivity of the sum rule to this quark-hadron duality ansatz. 

The lowest hadronic state is the kaon. Using the standard definition of the kaon decay 
constant 

(0\s-f^ 5 u\K + (q)) = iq„f K , (29) 
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one obtains the relevant hadronic matrix element of the pseudoscalar current: 

(0| ji 9) \K+(q)) = f K m\ , (30) 

so that the kaon contribution to the hadronic spectral density reads: 

Pk\s) = fK rn K^( m K - s). (31) 

The two heavier pseudoscalar resonances [35] are K\ = K (1460) and K 2 = /l(1830) with 
the masses = 1460 MeV and = 1830 MeV and total widths T Kl = 260 MeV and 
T k 2 — 250 MeV, respectively. These resonances are not yet well established, in particular, 
no experimental errors are attributed to their masses and widths. In any case, it seems 
plausible that the hadronic spectral density in the pseudoscalar channel with strangeness is 
dominated by the kaon and K X2 resonances, making this channel less complicated than the 
scalar channel where the strong Kn scattering in S-wave ( J p = 0 + ) demands a dedicated 
analysis (see e.g., [17,36]). 

A detailed analysis of the hadronic part in the pseudoscalar sum rules (in both FESR 
and Borel versions) is presented in [18], employing the hadronic spectral density where the 
contributions of two resonances K X2 with finite widths are simply added to the ground- 
state term of the kaon. Here we adopt the same ansatz for the hadronic spectral density 3 
in the sum rule (6): 



fKi m Ki B Ki( s ) » 

£= 1,2 


(32) 


where BkA s ) are the finite-width (Breit-Wigner type) replacements of the ^-function in 
the spectral density for K i j2 : 


S(m 2 K . - s) 


BkAs ) 


}_ f _ r Ki m Ki _ 

vr \(s- m 2 K ) 2 + (r Ki m K A 2 


(33) 


In [18] using FESR, the decay constants fxi and fx 2 of K x and K 2 resonances (defined 
similarly to fx) were fitted. As anticipated from ChPT, small values, in the ballpark of 
20-30 MeV for both f Kl and fx 2 were obtained. On the other hand, clue to the large mass 
multiplying these constants, the effects of K x and K 2 are quite noticeable in the hadronic 
part of the sum rule, hence, one has to avoid too large values of M 2 . We will use the 
estimates of fK 1 ,i< 2 from [18] as hadronic inputs in our numerical analysis of (6). 

Further improvements of the hadronic ansatz are possible, but they our beyond our 
scope here. In particular, it seems important to investigate the role of multiparticle states 
in the hadronic spectral function, starting from the two-particle states K*ir,Kp. In [18] 
it is assumed that multiparticle effects are at least partially taken into account in the 
finite widths of K X 2 . One usually neglects the possible contributions of the nondiag¬ 
onal transitions to the hadronic spectral function, e.g., intermediate states of the type 

3 For a different hadronic ansatz including I\*tt state explicitly, see [37]. 
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(0\j^\K)(K\K*Tr)(K*n\Ki)(Ki\j^\0). The analysis of the light-quark vector channel 
( J 1 ’ = 1 _ ) without and with strangeness (see e.g.,[38]) indicates that the effects of mix¬ 
ing between separate resonances via intermediate multiparticle states could be noticeable. 
Here, adopting the ansatz (32) we tacitly assume that the total widths of Ad ,2 account for 
the dominant contributions of multiparticle states. In order to estimate the influence of 
this effect, we will also consider a version of the hadronic spectral density (32) with the 
total widths of Ab j2 set to zero, interpreting the difference of the result with and without 
the widths as a rough estimate of the uncertainty due to multiparticle hadronic states. 


5 Inputs and numerical results 


For the running of the strong coupling a s and of the quark masses in MS scheme we use 
the four-loop approximation and employ the numerical program RunDec described in [39]. 
The reference value for the quark-gluon coupling is taken as a s (mz ) = 0.1187 [35]. The 
alternative choice a s (m T ) = 0.334 [35] produces a small difference which we include into 
the overall counting of uncertainties. We do not attempt to fit the u- and d-quark masses 
from the analogous sum rules, and simply take the current (non-lattice) intervals from [35]: 
m u (2 GeV) = (1.5 - 5.0) MeV, m d {2 GeV) = (5.0 - 9.0) MeV. 

The renormalization scale in our numerical calculation is taken as ji = M, reflecting the 
average virtuality of perturbative quarks and gluons in the correlator. In order to study 
the scale dependence we also vary the scale within M 2 / 2 < /i 2 < 2 M 2 . The window of 
Borel parameter is taken as in [18], 2 < M 2 < 3 GeV 2 . This choice allows to avoid large 
nonperturbative effects, simultaneously keeping the excited state contributions reasonably 
small. 

The remaining input parameters used for the OPE of the correlation function are: the 
quark condensate densities taken from GMOR relation (uu) = — / 2 m 2 /(2(m u + m,i)) where 
f n = 130.7 MeV [35]; the ratio of strange and nonstrange condensates (ss) / (uu) = 
0.8 ± 0.3; the gluon condensate density (a s /itGG) = (O.O^lg^) 1 ®) GeV 4 . Finally, the 
dimensionful parameter for the quark-gluon condensate density is taken as ml = 0.8 ± 0.2 
GeV 2 (for a recent comprehensive review of condensates see [40]). 

First of all, we address the main question which interests us here, namely, how good 
is the convergence of OPE for the correlation function in N 4 LO, and how large is the 
numerical impact of the new 0(a 4 ) correction. For that we define the ratios: 


, r O=0,2) 

' n 


(. M 2 ) 


{[n( 5 )"(M 2 )]g-°’ 2) }°(° 


[n (5 )"(m 2 )]S 2) + [n( 5 )"(M 2 )]g;J 6) 


(34) 


for n = 0,1,2, 3,4, where the numerator contains the contribution of 0(cc”) to the Borel 
transformed correlation function. The analogous ratio for the nonperturbative contri¬ 
butions is r^ =4,6 ^(M 2 ), where the numerator contains only the power suppressed term 
[n (5) "(M 2 )] op£ 6 \ Altogether, 




M 2 [ GeV 2 ] 


Figure 1: Relative contributions to OPE of the correlation function n ( ^" (M 2 ) defined in 
(34), plotted as functions of the Borel parameter squared. The solid lines from up to down 
correspond to rn~ 0,2 ' > with n = 0,1, 2, 3, respectively, the dashed line to r^ -0,2 '* and the 
dotted line to r^ =4,6 ^. 


^2 r^ =0 ’ 2) 

n=0,l,2,3,4 


_)_ r (d=4,6) 


1 . 


Note that the dominant rri s -dependence in hb 5 )" (M 2 ) is due to the overall factor (m s +m u ) 2 
and largely cancels in r < f l ~ 0 ’ 2 ' > and r^ =4,6 \ (We use m s (2GeV) = 105 MeV in the suppressed 
terms for this numerical illustration). In Fig. 1 the ratios r G -0,2 ) and ?4 d=4 ’ 6 ) are plotted as a 
function of Borel parameter squared. The convergence is excellent, even beyond the region 
of the Borel parameter chosen for the sum rule analysis. In the central point M 2 = 2.5 
GeV 2 we obtain ri d=0,2) (2.5 GeV 2 ) = 52.4%, 28.3%, 14.4%, 4.0%, -0.3% for n = 0,1, 2, 3,4, 
respectively and rG =4 > 6 )(2.5 GeV 2 ) = 1.2%. We conclude that the currently achieved 
accuracy of the correlation function at large M 2 is quite sufficient for the applications, 
such as the quark mass determination. 

We then turn to the sum rule (6). The input parameters for the kaon contribution to the 
hadronic part are: fx = 159.8 MeV, nix = 493.7 MeV [35]. The masses and total widths 
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of K] and K 2 resonances [35] were already quoted in the previous section. Their decay 
constants are taken from [18]: f Kl /V 2 = (22.9 ± 2.4) MeV, _/% 2 /-v/ 2 = (14.5 ± 1.5) MeV, 
where the factor 1/a/2 accounts for the difference between the normalizations. Note that 
for consistency, we take the version of [18] obtained without ILM correction; furthermore, 
we added the uncertainties of f Kl , f K2 given in [18] in quadrature. The duality threshold 
adopted in our calculation is So = 4.5 ± 0.5 GeV 2 . The central value provides the best 
stability of the sum rule in the Borcl parameter interval 2-3 GeV 2 , whereas the spread of 
±0.5 GeV 2 is added to allow for some additional variation of the hadronic input. For the 
middle values of all parameters specified above, we calculate m s from (6) and obtain the 
central value presented in (1). 

The influence of the new 0(a 4 ) correction on the sum rule is estimated by repeating the 
calculation with the same input, but with the perturbative corrections up to O(cc^). The 
result for the central value turns out to be m s {2 GeV) = 107 MeV, only 2 MeV larger than 
in (1). We also checked the quality of the OPE in the sum rule. Isolating the OPE part in 
(6), that is, subtracting from l.h.s. of (6) the integral over Pq PE {s) on r.h.s., we calculated 
the ratios fn~ 0,2 \M 2 , s 0 ) and r (cZ=4,6 ^(M 2 , s 0 ) defined analogous to Eq. (34), where instead 
of [II (M 2 )]qpe the contributions to the subtracted correlation function 

/»oo 

M 4 [n ( 5 ) "(M 2 )] OPE - / dse- s / M2 p^ PE (s) 

J so 

are substituted. We obtain fn~°’ 2 \ 2.5 GeV 2 ,4.5 GeV 2 )=39.2%, 26.1%, 18.8%, 10.6%, 
3.7% for n = 0,1, 2, 3,4, respectively, and f( d=4 > 6 )(2.5 GeV 2 ,4.5 GeV 2 ) = 1.6% revealing 
again a good convergence. 

Being numerically small, the 0(a 4 ) correction is nevertheless important for achieving a 
better stability with respect to the variation of the renormalization scale /i entering the sum 
rule. To demonstrate that, we have calculated rn s {2 GeV) from the sum rules with O(a^) 
and 0(a 4 ) accuracy, varying p 2 /M 2 from 0.5 to 2.0. The results plotted in Fig. 2 clearly 
demonstrate the role of the new 0(a 4 ) correction in stabilizing the scale-dependence. 

To investigate separate theoretical uncertainties of the sum rule (6) in more detail, we 
group them into two categories: 

a) uncertainties related to the input parameters in the correlation function and in 
the sum rule: renormalization scale, difference between using a s (m z ) and a s (m T ), Borcl 
parameter, u- and d-quark masses, condensate densities; 

b) uncertainties caused by the hadronic input: the decay constants fx 1 arid fx 2 , the 
effective threshold so and the effect of switching off the total widths of A% 2 - 

Varying the input parameters in the QCD part of the sum rule within the limits spec¬ 
ified above, we End that the the largest uncertainty in the category (a) is caused by the 
scale variation (see Fig. 2), whereas the sensitivity to the Borel mass variation is less 
than ±0.1 MeV, and the dependence on the values of condensate densities is negligible. 
Adding separate uncertainties grouped in this category in quadrature, we obtain the inter¬ 
val (±6 MeV) \param, included in (1). 

To investigate the hadronic uncertainties grouped above in the category (b), the decay 
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Figure 2: Strange quark mass at the scale 2 GeV, calculated from the sum rule (6) as a 
function of the renormalization scale p 2 in the correlation function, varying /i 2 /M 2 from 
0.5 to 2.0 at M 2 = 2.5 GeV 2 . The solid (dashed) line represents the result obtained with 
0 (a).) (0(a))) accuracy. 

constants fx 1 , ftc 2 and the threshold so are varied one by one yielding ±5 MeV, ±3 MeV 
and ±3 MeV, respectively. To estimate the effect of multiparticle states in the sum rule the 
rn s calculation is repeated with the total widths of K X 2 hr (33) set to zero. The result for 
rh s increases by approximately 3 MeV, which we conservatively interprete as an additional 
uncertainty ±3 MeV. All these individual uncertainties are again added in quadrature to 
produce (±7 MeV)| /iadr in (1). 

The nonperturbative effects beyond OPE cannot be estimated without the knowledge of 
the instanton effects, which are beyond our scope and are therefore absent in (1). According 
to the estimate [18] one has to add ±9 MeV to the total budget of uncertainties. In 
Table 1 we compare our prediction with other determinations of m s . The values of the 
s quark mass obtained from the correlation functions (r decays, Borel sum rules and 
FESR ) are consistent with each other and with the lattice QCD results within still large 
uncertainties, the lattice results with rif = 3 being systematically lower. Our estimate (1) 
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Method 

m s { 2 GeV) 

Ref. 


[MeV] 


Pseudoscalar Borel sum rule 

105 ± 6 ± 7 

This work 


100 ±6 

[18] (no ILM) 

Pseudoscalar FESR 

100 ± 12 

[18] 

Scalar Borel sum rule 

99 ± 16 

[17] 

Vector FESR 

139 ±31 

[19] 


81 ±22 

[10] 

Hadronic r decays 

9«+ 5 + 16 

yo -3-18 

[11] 


104 ± 28 

[22] 

r decays © sum rules 

99 ±28 

[22] 


97 ±22 

[2] 

Lattice QCD (rif = 2) 

100 -130 

[4] 


101 ± 8 Ijf 

[5] 


76 ± 3 ± 7 

[1] 

Lattice QCD (rif = 3) 

86.7 ±5.9 

[3] 


87 ± 4 ± 4 

[6] 

PDG04 average 

80 -130 

[35] 


Table 1: Our estimate of m s (2GeV) compared with some recent determinations obtained 
with different methods. The error/uncertainty identification in the results taken from the 
literature can be found in the corresponding papers. 


is also consistent with the s-quark mass bound in 0(a*) obtained in [23]. 

6 Conclusion 

We have included the new O(a^) correction in the correlation function of the pseudoscalar 
strangeness-changing quark currents and calculated the s-quark mass from the resulting 
Borel sum rule. In future the same analysis should be repeated for the scalar Borel sum 
rule, for both pseudoscalar and scalar FESR and for the sum rules with nonstrange light- 
quark currents, yielding the u and d quark masses. 

Our main intention here was to investigate the role of the Ofaff) terms in the OPE 
and in the sum rule. We have found that the new correction is comfortably small, making 
OPE in this channel very reliable. Simultaneously, the addition of the O(af) contributions 
noticeably decreases the renormalization scale-dependence of the resulting sum rule. 

The QCD sum rules obtained on the basis of OPE still have a considerable room of 
improvement. While the nonperturbative effects beyond OPE can be kept under control by 
choosing the virtuality (Borel parameter) scale large enough, and using ILM-type estimates, 
there is still a lack of experimental information concerning the masses, total and partial 


12 




widths of the excited kaon resonances. The resonance K\ can be observed in r —» Kmrv T 
decays, and both K\ and K 2 probably also in hadronic B decays where the currently 
available statistics allows to isolate many light-quark resonances in the final states. With 
this information one would be able to built a hadronic spectral function in the pseudoscalar 
channel which is less dependent on duality ansatz, so that the accuracy of the hadronic 
part of the sum rule eventually becomes closer to the high precision achieved in the QCD 
part. 
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